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ABSTRACT 

We make a numerical study of gauge and Yukawa unification in super- 
symmetric grand unified models and examine the quantitative implications 
of fermion mass ansatze at the grand unified scale. Integrating the renor- 
malization group equations with ai(M z ) and a 2 (M^) as inputs, we find 
a 3 (M z ) ~ 0.111(0.122) for M SUSY = m t and a 3 (M z ) ~ 0.106(0.116) for 
Msusy = 1 TeV at one-loop (two- loop) order. Including b and r Yukawa 
couplings in the evolution, we find an upper limit m t ^200 GeV from Yukawa 
unification. For given 777^175 GeV, there are two solutions for /3, one with 
tan/? > m t /m b , and one with sin/3 ~ 0.78(m t /150 GeV). Taking a popular 
ansatz for the mass matrices at the unified scale, we obtain a lower limit on 
the top quark mass of 771^150(115) GeV for a 3 (Mz) = 0.11(0.12) and an up- 
per limit on the supersymmetry parameter tan/?^50 if a 3 {M z ) = 0.11. The 
evolution of the quark mixing matrix elements is also evaluated. 



I. INTRODUCTION 



There is renewed interest in supersymmetric grand unified theories (GUTs) Q to explain 
gauge couplings, fermion masses and quark mixings &-H . Recent measurements of the gauge 



couplings at LEP and in other low energy experiments [10,11 ] are in reasonably good accord 
with expectations from minimal supersymmetric GUTs with the scale of supersymmetry 
(SUSY) of order 1 TeV or below ||. Supersymmetric GUTs are also consistent with the non- 
observation to date of proton decay ||12[| . In addition to the unification of gauge couplings 
|TB| |, the unification of Yukawa couplings has been considered to predict relations among 



quark masses [|i4Hl6|1- With equal 6-quark and r-lepton Yukawa couplings at the GUT scale, 
the mb/m T mass ratio is explained by SUSY GUTs Pjjlql. With specific ansatze for the 
GUT scale mass matrices (e.g. zero elements, mass hierarchy, relations of quark and lepton 
elements), other predictions have been obtained from quark masses and mixings that are 
consistent with measurements |f|,|6|]7l,|llj,|18|j . The consideration of fermion mass relationships 
has a long history [0,0 aric ^ includes single relations and mass matrices ( "textures" ) without 
evolution pT| , p2| , and single relations and mass matrices with evolution |23 . 



Our approach is to explore supersymmetric GUTs first with the most general assump- 
tions, and then proceed to add additional GUT unification constraints to obtain more pre- 
dictions at the electroweak scale. The renormalization group equations (RGEs) used here 
are for the supersymmetric GUTs [U,^l with the minimal particle content above the super- 



symmetry scale and the standard model RGEs |26[ below the supersymmetry scale. In §11 
we explore the running of the gauge couplings in the supersymmetric model at the two-loop 
level and compare the results to those obtained at the one-loop level. Rather than try to 
predict the scale of supersymmetry {M SUSY ) which may be sensitive to unknown and model 
dependent effects like particle thresholds at the GUT scale, we choose two values of M SUSY 
to illustrate the general trends that occur. We also investigate the effects of the Yukawa 
couplings on the gauge coupling running which enter at two loops [|17]] and have often been 
neglected in the past. In §111 we explore the unification of Yukawa coupling constants. First 
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we consider the one-loop analytic solutions which can be obtained by neglecting the bottom 
quark and tau Yukawa couplings A& and A T relative to At in the RGEs. This serves as a useful 
standard for comparison with the two-loop results for smaller values of tan/5 (<< mt/ra&), 
and many of the general features of the solutions to the RGEs are already present at this 
stage. We then investigate the two-loop RGE evolution of the Yukawa couplings including 
the effects of A;,, A T , and Af. Analytic solutions are not available for the two-loop evolution, 
so we integrate the RGEs numerically. In §IV we investigate relations between Cabibbo- 
Kobayashi-Maskawa (CKM) matrix elements and the ratios of quark masses. We investigate 
two popular ansatze @J7|,[1(J for Yukawa coupling matrices at the GUT scale. Both of these 
ansatze agree with all existing experimental data, and this agreement is preserved at the 
two-loop level. We also integrate the two-loop evolution equations for certain CKM matrix 
elements and quark mass ratios in §IV. The two loop RGEs for both the minimal supersym- 
metric model and the standard model are given in the appendix. 



A consistent treatment to two loops in the running of the gauge couplings involves the 
gauge couplings gt and the largest Yukawa couplings A t , A;, and A T . From general expressions 



II. GAUGE COUPLING UNIFICATION 



24] , |25|| that are summarized in the appendix, we obtain the evolution equations 




(1) 
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( E O + 4 2 /2) ^ + s 2 s 2 2 + \g\gl + Sgjgt 



dX T X T 



16vr 2 



+^h 2 i + tf> (J*?? + 6^ 2 + 16^ 2 ) + l^gf 

- {22Xt + 5A 2 A 2 + 3A 2 A 2 + 3A* + 5A*} 



-E« + 3A 2 + 4A^ 

+ I ^(E(^ + < 2 /2). 4 + ^ 2 2 

±A 2 (- 2 -g 2 + lQg 2 )+X 2 T (jg 2 + Qg 2 2 
-{3X 2 X 2 b + 9Xt + 9X 2 b X 2 T + WXt} 



(3) 



(4) 



The various coefficients in the above expressions are also given in the appendix. The variable 
is t = \ii(h/Mg) where [i is the running mass scale and Mq is the GUT unification mass. 
The renormalization group equations of dimensionless parameters like the gauge couplings 
and Yukawa couplings are independent of the dimensionful soft-supersymmetry breaking 
parameters. 

We begin with the recent values of a em and sin 2 9 W at scale M z = 91.17 GeV given in 



the 1992 Particle Data Book [0.127 



127.9 ±0.2 



sin 2 6 W = 0.2326 ± 0.0008 



(5a) 
(5b) 



where 9 W refers to the weak angle in the modified minimal subtraction MS scheme [£8 
These values correspond to electroweak gauge couplings of 



a^Mz)- 1 = 58.89 ±0.11 
a 2 (M z )- 1 = 29.75 ±0.11 



(6a) 
(6b) 



For simplicity we initially set the supersymmetric scale M SUSY equal to the top quark mass 
m t and set all Yukawa contributions in Eq. (|IJ) to zero. Then evolving ot\ and a.2 from scale 
Mz up to scale m t , we have 



53 

a^mtY 1 = ol~l(Mz) + —ln{M z /m t ) 

OVTT 

a 2 (m t )- 1 = a 2 {M z y 1 - — hx{M z /mt) , 

07T 



(7a) 
(7b) 



We use the value M z = 91.17 GeV, neglecting its experimental uncertainty. 

Next, for a grid of a G and M G values, we evolve from the GUT scale down to the chosen 
m t scale and retain those GUT scale inputs for which Eqs. (6) and (7) are satisfied. We use 
the two-loop SUSY GUT unification condition a G = a\(M G ) = ct2(M G ). For the acceptable 
GUT inputs we also evolve the strong coupling a^(M G ) = a G down to scale m< and then use 
3-loop QCD to further evolve it to scale M z . The three-loop expression 

-i 



a 3 (/x) 1 



'2 ln N VH ^ ~ 



In In 



e 

A 2 



bob 



Q u 2 



b\ 



1 



In 



EL 

A 2 



with the bi given in Ref. [29"|, is iteratively solved to find A from a 3 (m t ). Eq. (||) is then 
evaluated for fi = M z to obtain a^(M z ). The resulting values for A for two representative 
values of q>z(M z ) are given in Table 1. 
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a 3 {M z ) 


A (5) 


A (4) 


A (3) 


0.11 


129.1 


188.3 


225.0 


0.12 


233.4 


320.2 


360.0 



Table 1: The QCD parameter A^) in MeV, 
where rif is the number of active flavors. 



We also investigate the effects of taking a supersymmetry scale higher than mj. Be- 
low M SUSY , the RGE are similar to the non-supersymmetric standard model. A linear 
combination of the Higgs doublets is integrated out of the theory at M SUSY leaving the 
orthogonal combination $(sm) = ®d cos P + ^« s * n ft coupled to the fermions in a way that 



depends on tan/5 P J30|J31[] ; this combination results from the assumption that the three 
soft-supersymmetry breaking parameters in the Higgs potential can be equated to M SUSY . 



We use the two- loop RGEs |26| for the standard model, matching the couplings at M SUSY . 
Taking a single SUSY scale is an idealized situation since in general the supersymmetric 
particle spectrum is spread over a range of masses ||. Without further assumptions we 
cannot predict this spectrum. Given that such uncertainties exist, the predicted range for 
«3 should be taken to be representative only. 

The ranges of a^ 1 and M G parameters obtained from the procedure outlined above are 
presented in Fig. 1 for one-loop and two-loop evolution with the choices M SUSY = m t and 
M SUSY = 1 TeV. The shaded regions denote the allowed GUT parameter space. The two- 
loop values obtained for a G and M G are higher than the one-loop values and consequently 
a^iMz) is higher for the two- loop evolution. Note that raising the SUSY scale from m t to 1 
TeV lowers M G and a G ; hence a^{Mz) decreases as well. 

Figure 2 shows the corresponding results of the two-loop evolution over the full range of 
fi. We find the ranges for a^(Mz) with m t = 150 GeV shown in Table 2. 
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Msusy 


one-loop 


two-loop 


m t = 150 GeV 


0.1112 ±0.0024 


0.1224 ±0.0033 


1 TeV 


0.1065 ±0.0024 


0.1161 ±0.0028 



Table 2: Ranges obtained for as(Mz) from the 
input values a em and sm 2 9 w . 



The two-loop values of ots(Mz) are about 10% larger than the one-loop values. The effect of 
the higher SUSY scale is to lower a^{Mz) by about 5%. 

Inclusion of Yukawa couplings in the two-loop evolution also lowers the value of a%{Mz) 
somewhat. For example setting At = A& = A T = 1 at the GUT scale, we obtain a two- loop 
value of a 3 {M z ) = 0.1189 ± 0.0031 for M SUSY = m t . 

The effects on the gauge couplings of including the Yukawa couplings in the evolution 
are rather small for Yukawa couplings in the perturbative regime, justifying their neglect in 
most previous analyses; for large values of tan (3 the changes in the gauge couplings due to 
inclusion of Yukawa couplings can be a few percent. 

The experimental situation regarding the determination of a 3 is presently somewhat 
clouded flOf , with deep inelastic scattering determinations in the range of the one-loop 
calculations in Table 2 and LEP determinations similar to the two-loop results of Table 2. 

There are other uncertainties not taken into account here, due to threshold corrections 
from the unknown particle content at the heavy scale [|32} -|34|1 , which can also change the a 3 
values obtained above. These corrections are mo del- dependent so we have not attempted 
to include such contributions. However recent analysis have shown that the constraints 
from non-observation of proton decay greatly reduce the potential uncertainties from GUT 



thresholds ||17|j35| . 
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III. YUKAWA UNIFICATION 



A. One-loop analytic results 

The unification of Yukawa couplings first introduced by Chanowitz, Ellis and Gaillard |I3 
has been reconsidered recently P|]fi|,|7|,|r7|,ftm . The GUT scale condition X b (M G ) 



r7|,g. The GUT scale condition X b (M G ) = X T (M G ) 
leads to a successful prediction for the mass ratio mb/m T provided that a low energy super- 
symmetry exists H. The b to r mass ratio is given by 



m b _ r] b 



R b / T (m t ) 



(9) 



where 



Rb/ T (m t ) = 



X b (m t ) m b {m t ) 



X T {m t ) m T (m t ) 
is the b to r ratio of running masses at scale mt and 



(10) 



Vf 



Vf 



m f( m f) 
m f (m t ) 



m/(lGeV) 



if rrif > lGeV 



if rrif < lGeV 



(12) 



mf{m t ) 

is a scaling factor including both QCD and QED effects in the running mass below m t . We 
have determined the T]f scaling factors to three- loop order in QCD and one- loop order in 
QED. The QCD running of the quark mass is described by 

7(A 



mM = m q (2b a 3 y o/bo 



+ - 



7i 
bo 



lob 



1 + 

\ 2 



7i 
bo 



h 2 



Oi 3 



bl 



Uo &g bl , 



«3 + 



(13) 



where the anomalous dimensions 70, 71 and 72 are given in Ref. [36]. The scale-invariant 
mass m q cancels in the ratio in Eq. (|Tl|) . The one- loop QED running from scale // to scale 
/x introduces modifications 

7o QED /"o QBI5 



m/(/i) = m f (p') 



a{jj) 



(14) 



where the QED beta function and anomalous dimension are given by [[33 



l8 ED = -3Q} , (16) 

and the sums run over the active fermions at the relevant scale. The dependence of the QCD- 
QED scaling factors t] on a^iMz) is shown in Figure 3; these factors increase as a^iMz) 
increases. 

We note that the physical top mass is related to the running mass by WA 
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rrit hys = m t (m t ) 



-a 3 (m t ) + 0(a 2 3 ) 
on 



(17) 



The effects of the top quark Yukawa At can be studied semi-analytically at one-loop 
neglecting the effects of the bottom and tau Yukawa couplings A& and A T in Eqs. (|I|) and 
(Q), which is a valid approximation for small to moderate tan/3 (i.e. tan/3^10). Following 



Ref. M we find [38 



m b V 1/2 Vb , 1fi , 
— = y , 18 

m T x r] T 



where x(fi), T](fi) defined by 



x(») = (a G /«i(/ i )) 1/ V G /a 2 (rf) 3/2 , (19) 

r M a } 

y(^) = exp|-— J A t 2 (//)cnn/i' , (20) 
Vfr) = II (ao/oiCA*))^ , (21) 

i=l,2,3 



are to be evaluated at fx — m t in Eq. fll8P. Henceforth x, y, rj shall be understood as being 
evaluated at scale mt when an argument is not explicitly specified. Typical values of these 
quantities obtained in Ref. [18| are x = 1.52, y = 0.75 — 0.81, 77 = 10.3 for a bottom mass 
given by the Gasser-Leutwyler (GL) QCD sum rule determination mj = 4.25 ±0.1 GeV [B9| 



taken within its 90% confidence range and = 0.111. The quantity y gives the scaling from 
Mq to m t that arises from a heavy quark, beyond the scaling due to the gauge couplings. 
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The factor y(m t ) is constrained to lie in a narrow range of values by Eq. flTq). The integral 
in Eq. ( f20"D is crucial in explaining the mb/m T ratio. In fact if A f is neglected then y = 1 and 
the rrib/m T ratio is found to be too large. 

For a given value of mt, there exist two solutions for tan/?. This fact can be understood 
qualitatively by studying the one-loop RGE for R h / T = A&/A T . 

^ = ^(-E*ft 2 + A? + 3^-3Ay . (22) 

For small tan {3 the bottom and tau Yukawas do not play a significant role in the RGE, and 
any particular value for mb/m T is obtained for a unique value of \ t (m t ), which corresponds 
to a linear relationship between m t and sin j3, 



m t v 



sin/3 v/2 
where v = 246 GeV and 



1/2 



\ t (m t )=7rvi 1-y 12 , (23) 



77(//)dln// . (24) 



The numerical value for / from Ref. |18| is 113.79 for m t = 170 GeV. For large tan/3, where 
the effects of A;, and A r on the running Yukawa couplings can be substantial, an increase 
in A;, can be compensated in the RGE by a decrease in X t . Hence, for increasing tan/3, the 
correct prediction for mb/m T is obtained for decreased values of the top quark Yukawa. Thus 
there is a second solution to the RGE for R\,It with a large value of tan/3. The inclusion of 
the two-loop effects does not alter these observations. 
The one-loop RGE for R s /^ = A s /A^ 

^ = {!£(-£«<*?)■ (25) 

is similar to Eq. (E2), except that it receives no contribution from the dominant Yukawa 
couplings At, Af,, and A T . When the value R s /^(Mg) = 1/3 is assumed at the GUT scale, the 
prediction at the electroweak scale is 

m s 1 77 1 / 2 r] s 



m M 3 x rjf, 
10 



(26) 



Notice that this equation does not include the scaling parameter y because the top quark 
Yukawa does not affect the running of the second generation quarks and leptons. This 
relation for m s /m M is in good agreement with the experimental values, but it is not as 
stringent as the mb/m T relation due to the sizable uncertainty in the strange quark mass. The 
result m s /m M = 1.54 was obtained in Ref. []!£[] , to be compared with the GL determination 



m s /m^ 



1.66 ±0.52. 



A popular strategy is to relate the mixing angles in the CKM matrix to ratios of quark 



masses, taking into account the evolution from the GUT scale in non-SUSY J4C§ or SUSY 



| ill models. For example, one popular GUT scale ansatz is \V cb \ 
a GUT boundary condition on R c / t = A c / At of 



y/R c/t (M G ) = \V cb {M G ) 



m c /m t which requires 



(27) 



The one-loop SUSY RGE for R c / t is 

dR c /t _ 
dt ' 



R 



c/t 



16tt 2 



3A? + A^ 



(28) 



The corresponding one- loop SUSY RGE for the running CKM matrix element \V cb \ is ||T 

d\V cb \ 



dt 



\Vcb\ 

16tt 2 



K + K 



(29) 



The pure gauge coupling parts of the RGEs are not present in Eqs. ( p8|) and (p9|) since R c / t 
and V cb are ratios of elements from the up quark Yukawa matrix and the down quark Yukawa 
matrix. 

Neglecting the non-leading effects of X b , the one-loop results of Ref. || at the electroweak 
scale obtained from evolution are 



\V cb (m t )\ = iV^M^ly- 1 R c/t (m t ) = R c/t (M G )y- 



(30) 



or equivalently using Eq. ( p7| ) 



\V cb (m t )\ 



lym c 
r] c m t 



(31) 
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Since y is already well constrained by the 6-mass relation of Eq. (|18D (for the one-loop 
value of a 3 (M z ) = 0.111), Eq. (|31|) requires that m t must be large in order that |V^,| falls 
in the experimentally allowed range 0.032 — 0.054 (and even then \ V c b\ is found to be at the 
upper limit of its allowed range). If, however, we use a larger value of [Mz) indicated 
by the two-loop equations, say 0.12, then r\ c increases by about 14%, as shown in Figure 
3. Furthermore the increased values of the scaling parameters r\ and rjb require about a 9% 
decrease in y to explain the mj/m T ratio in Eq. ([Tj|). The resulting \V c b\ is reduced by about 
12% and is then closer to its central experimental value. Of course, a consistent treatment 
at the two-loop level requires the two-loop generalization of Eq. ( |31~D obtained by solving the 
full set of RGEs. One of the questions we will address subsequently is for what values of m t 
and tan/3 can the \V c b\ and the mf,/m T constraints be realized simultaneously. 

The predictions above are all based upon the assumption that the couplings remain in 
the perturbative regime during the evolution from the GUT scale down to the electroweak 
scale. Otherwise it is not valid to use the RGEs which are calculated order by order in 
perturbation theory. One can impose this perturbative unification condition as a constraint. 
For rrib at the lower end of the GL QCD sum rule range 4.1 — 4.4 GeV the top quark Yukawa 
coupling at the GUT scale, X t (M G ), becomes large, as can be demonstrated from analytic 
solutions to the one-loop RGEs in the approximation that A& and A r are neglected compared 
to At (valid for small to moderate tan/3). 

The top quark Yukawa at the GUT scale is given by 



4vr 2 

A,(M G ) 2 = i" 



4-1 

yll 



(32) 



Taking |8f a 3 (M z ) = 0.111 and m b = 4.25 GeV and m t = 170 GeV gives X t (M G ) = 1.5. 



Larger values of a 3 (M z ) lead to increased r\ h via Eq. (|IID giving smaller y in Eq. (|1^) and 
a correspondingly larger value of X t (M G ). The quantity X t (M G ) is plotted versus a 3 (M z ) in 
Figure 4. Larger values of a^iMz) ~ 0.12 can yield Xt(M G )^,3 that cast the perturbative 
unification in doubt. Keeping the gauge couplings fixed and varying mt,, one sees that smaller 
values of rrib also yield larger values of X t (M G ). 
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The scaling parameter y is manifestly less than one by Eq. (|20|) since > in the region 
m t < n < M G . This implies an upper limit on rn b in Eq. ( fl8| ) of 

^01^ , (33) 
m T x r] T 

B. Two-loop numerical results 

When the two-loop RGEs are considered, analytic solutions must be abandoned, but 
the same qualitative behavior is found in the numerical solutions. Furthermore, there is 
now the possibility that the bottom quark Yukawa coupling at the GUT scale becomes 
non-perturbative for large values of tan/?. In our analysis we solve the two-loop RGEs of 
Eqs. ([l|-|4]) numerically |12 1 , retaining all Yukawa couplings from the third generation. 



First we choose a value of a^iMz) that is consistent with experimental determinations 
and the preceding one-loop or two-loop evolution of the gauge couplings in the absence of 
Yukawa couplings. Specifically we take as(Mz) = 0.11 or a^{Mz) = 0.12, to bracket the 
indicated a^(Mz) range. For each particular a^iMz) we consider a range of values for tan/? 
and m b {m b ). For each choice of a^Mz), tan/3, m;, we choose an input value of m t . The 
Yukawa couplings at scale m t are then given by 

y/2m t (m t ) V2m b (m b ) V2m T (m T ) 

M m t) = ^-5— > M m t) = ir ' M™t) = 7T ' ( 34 ) 

v sin p r] b v cos p r] T v cos p 

and the ai(m t ) are determined by Eqs. (7) and (|8|) from the central values in Eq. (6) We take 



lT| m T = 1.784 GeV. The running of the vacuum expectation value v between the fermion 
mass scales is negligible for the range of fermion masses considered here |J. Starting at the 
scale m t , we integrate the RGEs to the GUT scale, defined to be the scale at which ai(^) 
and «2(/i) intersect. We then check to see if the equality A&(M G ) = X r (M G ) holds to within 
0.01%. If the b and r Yukawas satisfy this condition, the solution is accepted. If not, we 
choose another value of m t and repeat the integration. Since our primary motivation here 
is to study the influence of the a^iMz) value on the Yukawa couplings, we do not enforce 
the requirement that a^{M G ) is equal to a\(M G ) and 02 (M G ). Nevertheless the equality of 
«i, «2, and «3 at M G is typically violated by ^ 4% (2%) for as(Mz) = 0.11 (0.12). Such 
discrepancies could easily exist from threshold effects at the GUT scale []3 
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We also explore the effects of taking the SUSY scale above m t . We proceed as described 
above, integrating the following two-loop standard model RGEs numerically from the top 
mass to the SUSY scale: 



dgi 

dt 



9i 
16tt 2 



b? M 9* + 



16tt 2 



u'=i 



j=t,b,T 



(35) 



A* 



16tt 2 



dA h 



A?, 



16tt 2 



+ 



1 /1187 
16^2 I 



G()0 9l - "J#2 - !08^3 - 20^?^2 + ^£l#3 + 9 #2#3 



/223 2 135 o 2 \ l2 /43 2 9 2 

+^Y 4 (S) - 2A (3A 2 + A 2 ) 
+^ - \\\\\ + ^Xt 
+ Y 2 (S) Qa 2 - ?A 2 ) - XA (S) + ^A 2 



\ 2 



+ 



[-£« + ^A 2 -^A 2 + Y 2 (S)] 

i ( - - t 52 - 108534 - i^ 2 + §^ 2 + 9 ^ 

+^Y 4 (S)-2\ (A 2 + 3A 2 ) 

5.o.n 11 



187 
80 



15 
135 
~16 



<?! + 16(? 2 ) A 2 



+ 3 A 4 



+Y 2 (S) Qa 2 - ^A 2 ) - US) + V 



(36) 



(37) 



dA T 



A r 
16tt 2 



+ 



cf M S 2 + ^A 2 +Y 2 (S) 



1371 4 23 4 27 2 2 

#1 - "J#2 + 20^1^2 



16tt 2 I 200 
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387 2 135 2 \ x2 5 V L , 



80 



16 



6AAf 



+\K-\y2{s)\1-x,{s) + \\* 



(38) 



dA 
~dt 



16tt 2 



9/3 
4 V25 



Here 



at + \glal + at) - (J<? 2 + % 2 2 ) A + 4F 2 (S)A - ah(S) + m 2 } 

1 ( „ ,^ /3 o ,\ l2 / 73 4 117 2 2 1887 A , 



305 fi 867 2 4 1677 

H 99 9799 

8 y2 120 1 2 200 

-64^ 2 (A 4 + A*) 



4 2 



3411 
1000 



9i 



--g\{2\t - A 4 + 3A 4 T ] - -g 4 2 Y 2 (S) + 10AF 4 (5) 
+\ 9 l [(-^i 2 + 21^ 2 ) A 2 + {^- g l + 9, 2 ) A 2 



-24A 2 F 2 (,S) - \H{S) + 6AA 2 A 2 



+20 
-12 



3A t 6 + 3A° + \ % T 



\ ^b + ^i^i 



(39) 



Y 2 {S) = 3 A 2 + 3 A 2 + A 2 , 



(40) 



Y*(S) = \ [3 E <fW + 3 E <fW + E 



(41) 



X4(5) 



3A 4 + 3A 4 + A 4 - 



(42) 



H(S) = 3A 4 + 3A 4 + A 4 , 



(43) 



and the coefficients a SM , b SM and c SM are given in the appendix along with the full matrix 
structure. 
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The initial values for az{Mz), vn b and m t are chosen as before; in addition we are required 
to specify the initial value of the quartic Higgs coupling A at scale m t . The Yukawa couplings 
at scale m< are 

^2m t {m t ) y/2m b (m b ) V2m T (m T ) 

M m t) = , h{m t ) = , M m t) = > (44) 

v rj b v rj T v 

and the ai(m t ) are given by Eqs. (6)-(8). After integrating to the SUSY scale we require 
that the matching condition 

KMsusy) = \ (^9l{M+ USY ) + g 2 2 (M+ USY )j cos 2 2(3 , (45) 

is satisfied to within 0.1%. This condition [§],|31| results from integrating out the heavy Higgs 
doublet at M SUSY . Below this scale only a Standard Model Higgs remains with its quartic 
coupling given by Eq. (l45|). We also apply the matching conditions 



9i( M SUSY) 


= 9i( M SUSY) 




(46) 


HMsusy) 


= MMtusY) 


sin j3 , 


(47) 


MMsusy) 


= HM^usy) 


cos j3 , 


(48) 


K{Ms USY ) 


= K{Mg USYi 


cos/3 . 


(49) 



If Eq. ( f45|) is not satisfied we choose another input value of X(m t ) and repeat the process. We 
allow tan (3 to span a wide grid of values. After obtaining a satisfactory value of A that meets 
the boundary condition above, we integrate the two-loop SUSY RGEs to the GUT scale, 
defined by the equality «i(M G ) = a 2 (M G ). At the GUT scale we require X b (M G ) = A T (M G ) 
to within 0.1%. If this condition is not met, we repeat the entire process, choosing other 
initial values for m t and A. 

The parameter (3 also runs in going from the SUSY scale to the electroweak scale [ |3Tf . 
However this effect is small and we neglect it here. 

In Figure 5 the resulting contours of constant m b are given in the mt, tan (3 plane |]4"|,|TT 



for the choices of a^{Mz) = 0.11 and 0.12 and the supersymmetry scales M SUSY = m t and 1 
TeV. The contours shown are m b = 4.1,4.25,4.4 GeV (corresponding to the central value of 
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m& and its 90% confidence range from the GL QCD sum rule determination) and nib = 5.0 
GeV (representing a typical constituent 6-quark mass value). For a given m& and m^l75 
GeV, there is a high solution and a low solution for tan /3 as anticipated in §IIIa. Thus, once 
rrit is experimentally known and the choice of nib resolved by other considerations (such as 
the CKM matrix elements addressed subsequently), the assumption of Yukawa unification 
at the GUT scale will select two possible values for tan (3. For example for m t = 150 GeV 
and nib = 4.25 GeV, the solutions with a^(Mz) = 0.11 are 



tan/3 = 1.35 or tan/3 = 56 . 

For m t &l75 GeV the low solution is well-approximated by 

nit 



(50) 



sin/3 = 0.78 



150 GeV 



(51) 



Such knowledge of tan/? would greatly simplify SUSY Higgs analyses pi] . Without imposing 
any other constraints, the top quark mass nit can be arbitrarily small. 

The plots rise very steeply for the maximal value of m t . This results because the linear 
relation exhibited in Eq. Q2"5|) and in the plot in Ref. [fTBl between m t and sin (3 is mapped 



into a vertical line for sufficiently large tan/3 (^10). The deviation of these contours from 
being strictly vertical results from the contributions of \ and A T to the Yukawa coupling 
evolution. 

An upper limit on m t is determined entirely by the nib/m T ratio. We find the m t upper 
limits shown in Table 3 for the two choices of a^iMz)- It is interesting that the predicted 
upper limit for m t coincides with that allowed by electroweak radiative corrections fTT|j. 



Msusy 


a 3 (M z ) 


0.11 


0.12 


m t 


187 


193 


1 TeV 


192 


199 
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Table 3: Maximum values of m t (m 4 ) in GeV consistent with 
the 90% confidence levels of the m^rrib) values of GL. 



Our contours of rrib/m T in Fig. 5 have about a 10% higher than those given in Ref. |T7| 
presumably because they employed the one- loop QCD results for the scaling factors r/j with 
the two-loop expression for «3 rather than the three-loop QCD for both r]f and 0:3 that we 
use here. 

As a^{Mz) gets larger, smaller values oiy are needed to obtain obtain the correct mb/m T 
ratio. In turn larger values of A 4 (/i) are needed via Eq. fl2T|). For a 3 (M^)^0.12 and m&£s4.2 
GeV, the value of Af(/z) near the GUT scale can be driven into the nonperturbative regime. 
In Figure 6 we show the values of At(Mc) and \b(Mc) obtained for the solutions in Fig. 
5. Fixed points in the quark Yukawa couplings exist at A ~ 1, so a Yukawa coupling only 
slightly larger than the fixed point at the scale m t can diverge as it is evolved to the GUT 
scale. For large values of the Yukawa couplings the two-loop contributions to the RGEs 
contribute a fraction x of the one-loop contributions when 

A ( = ^!^6VS, (52) 

A t = ^|^«6.3V5, (53) 

as can be deduced from Eqs. (Q) and @. When x ~ 1 we are clearly in the nonperturbative 
regime. If we adopt the criteria that the two-loop effects always be less than a quarter of 
the one-loop effects, then X t and A& are nonperturbative when they remain below 3.3 and 
3.1 respectively all the way to the GUT scale. This is true for all of the curves presented in 
Figure 6, except for the rrib = 4.1 GeV contours for a^iMz) = 0.12; hence the exact position 
of this contour cannot be predicted with accuracy. 

In Figure 7 we show the evolution of the Yukawa couplings from the SUSY scale to the 
GUT scale. The nonperturbative regime for the case discussed above occurs only near the 
GUT scale. 
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In some SO(10) GUT models the top quark Yukawa coupling At is unified with the A& and 
A T at the GUT scale. Imposing this constraint selects a unique value for tan/? and m t . This 
solution is given by the intersection of X t (Mc) and A&(Mg), which occurs for large tan/3^50: 
see Fig. 6. 

One could also consider the unification of the Yukawa couplings at some scale other than 
that at which the gauge couplings unify [p||l7|. Since R b / T increases as it evolves from the 
GUT scale to the electroweak scale, Yukawa unification at a scale larger than the gauge 
coupling unification scale gives a larger mb/m T ratio. 

The authors of Ref. 0] predict the light physical Higgs mass rather precisely. However 
this prediction is related to their assumption (and the one we use here) that the heavy Higgs 
doublet is integrated out at M SUSY . This means that the heavy physical Higgs bosons have 
masses Mh ~ Ma ~ Mff± ~ M SUSY » Mz, which requires that the light Higgs mass is 
close to its upper limit. The relation of sin/5 to m t then fixes the one-loop corrections to the 
light Higgs mass. 



IV. FERMION MASS ANSATZ 

By assuming an ansatz for Yukawa matrices at the GUT scale and evolving these matrices 
down to the electroweak scale, predictions can be obtained for the quark and lepton masses 
and the CKM matrix elements 0,^,^,|l6j. Much work has been done on individual relations 
such as \Vud\ ~ \Jm s /md and \V c b\ ~ \Jm c /m t which are imposed at the GUT scale as 
described in §111. Recently interest has been revived in models that involve several such 
relations, leading to a number of predictions for quark masses and CKM matrix elements at 
the electroweak scale |4]|^-|8|]. The relations evolve according to RGEs, and the main effects 
are determined by the largest couplings. For moderate values of tan/3 (i.e. tan/3^10), 
these are the gauge couplings and the top quark Yukawa coupling X t . For large values of 
tan/3(~ m t jmb) the effects of A& and A r can also be significant. Various individual relations 
at the GUT scale such as \V c b\ ~ Jmc/mt can be satisfied for certain choices of these Yukawa 
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couplings. The remarkable aspect of these fermion mass ansatze is that many relations can 
be made to work at one time. We shall concentrate in this section on two predictive ways 
of generating mixing between the second and third generations which put those mixing 
contributions entirely in the up quark Yukawa matrix ||,|6|,[16| or entirely in the down quark 
Yukawa matrix ]7[|. 



20 



A. The HRR/DHR Model 

Harvey, Ramond and Reiss |I6| proposed that the Yukawa matrices at the GUT scale 
have the form 



(54) 
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(55) 



These matrices incorporate both Fritzsch zeros pi and the Georgi-Jarlskog relation PT 



between down quark and charged lepton matrix elements. This relative factor of three has 
been realized in Higgs models with certain vacuum breaking patterns. HRR obtained the 
above ansatz using a 10 and three 126 Higgs multiplets in an SO(10) GUT model to obtain 
various relationships between CKM matrix elements and quark masses. The GUT ansatz of 
Eqs. (??) and ([55"D is also the basis for the recent RGE analysis by Dimopoulos, Hall and 
Raby ||. Henceforth we shall refer to this ansatz as the HRR/DHR model. It yields the 
relation \V c b\ : 



'X c /Xt at the GUT scale. 
Renormalization group evolution generates non-zero entries in the above Yukawa matrices 
and also splits B\ = U23 and B 2 = U32 to give the matrices at the electroweak scale of the 
form 
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(56) 
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(57) 
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The quantities A, D and D' are equivalent to At, A& and A r respectively up to subleading 

corrections in the mass matrix diagonalization. The one-loop solutions || to leading order 

in the hierarchy can be obtained analytically neglecting A& and A T . The one- loop results for 

the CKM elements at the scale mt are 

rj s m d rj c m u 



\V US 
\V cb 

v ub 



+ 

rj d m s r) u m c 



m s r] c m u m d 

+ 2a / cos ( 

V VdVum s m c 



1/2 



V, 



cb 



iym c 

rj c m t 
IVcm u 



(58) 
(59) 
(60) 



where r]i(m t ) is defined by Eq. ( [TT]) and y(m t ) by Eq. (|20|). The angle is a priori arbitrary. 
The down-type quark masses are related to the corresponding lepton masses by 

(61) 
(62) 
(63) 



m d 



m b = y- 



V 1/2 Vdn 
3m e 

X 7] e 

x Vf, 3 

i ii 
V 1 Vb 



X T) T 



Using the general expressions for the two-loop RGEs given in the appendix and keeping 
only terms unsuppressed by the hierarchy, one obtains Eqs. (p])-(||) as well as 
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(64) 
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(67) 



Notice that since l/B 2 dB 2 /dt= l/X t dX t /dt, the ratio B 2 /X t is constant over all scales and is 
in particular equal to its value at the GUT scale (B 2G / X t (M G )) . 

With these RGEs we can include the additional experimental constraints from the charm 
mass m c and the CKM matrix element \V c b\ to determine the allowed region of the HRR/DHR 
model in the m*, tan/3 plane. An analysis at the one-loop level neglecting X b and A r relative 
to X t was presented in Ref. [ITS . 



The Yukawa matrices are diagonalized by unitary matrices V£, V^, V^, Vj* so that 
Udiag = vfXJV** and D dia s = V d L ~DV d R] . The CKM matrix is then given by V CKM = 
V£V d \ We define a "running" CKM matrix by diagonalizing the Yukawa matrices U and 
D at any scale t. We find that A c /A t and \V c b\ are described in terms of the Yukawa matrices 
by 
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A, 
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(68) 



Vh \ = ?±- 6 -* 



Xt X b 



(69) 
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with 



— = rj c R c / t (m t ) ■ 
m t 



(70) 



To leading order in the mass hierarchy, the ratio R c /t is given by the ratio of eigenvalues 
of the 2x2 submatrix of U in the second and third generations while V c b is given by the 
difference in the rotation angles needed to rotate away the upper right hand entry in the 
submatrices of U and D. Given that the mass hierarchies exist, there is a simple iterative 
numerical procedure for diagonalizing the mass matrices U and D and obtaining the CKM 
matrix. We have checked that the corrections to the above formulas from contributions 
subleading in the mass hierarchy are small. 

It is straightforward to derive the resulting renormalization group equations from 



Eqs. 

dR c /t _ R c /t 



dt 



16tt 2 



d\V c 



cb\ 



w, 



cb\ 



dt 



16tt 2 



(3A t 2 + A| 

+ ^ (A t 2 (~<7 2 + Ggl) + \\ 2 b gl - (l3A{ + 2A?Ag + 5A* + A 2 A; 
(A? + A^ 



+ 



-Xhl + - (5A^ + 5A* + 



167r 2 V5 l " J1 5 

The corresponding evolution equations in the Standard Model are given by 
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yA< + ^Xj + H A £ + lx 2 X + \xlXl) j , (74) 
The evolution equations in Eqs. (|73|)-(|74"|) are obtained from the two-loop RGEs of the 



standard model given in Ref |2B[] and in the appendix. 

In the supersymmetric model \V c b\ increases with the running from the GUT scale to the 
electroweak scale |)41|| ; this is evident at the two- loop level in Eq. ([72]) . The opposite behavior 
occurs in Eq. ( [74]) for the nonsupersymmetric Standard Model where \V c b\ decreases as the 
running mass decreases [4(J. Fig. 8 shows the running of \V C b\ for the cases M SUSY = m t 



and 1 TeV. In contrast to \V c b\ the ratio R c u increases monotonically as the running mass 
decreases in both the Standard Model and supersymmetric model cases. 

We stress that Eqs. (|71|) and (ff4]) are the correct evolution equations regardless of the 
fermion mass ansatz at the GUT scale. Changing the ansatz just changes the boundary 
conditions at the GUT scale (terms subleading in the mass hierarchy differ between models, 



but this is a negligible effect). In a model for which the relationship \ V c b\ = J X c / X t holds (as 



in the HRR/DHR model), this boundary condition is y-R c / t (Mc) = |K;&(Mg)|. In Giudice's 
model, to be described below, the mixing between the second and third generations arises 
in the down quark Yukawa matrix alone, and so in his model R c / t and \V c b\ are unrelated at 
the GUT scale. 

In our analysis of the CKM constraints we proceed as in the discussion of the calculation 
for Figure 5. We numerically solve the two-loop RGEs as given by Eqs. ([l|)-(|4]), (|7l]) -(|72]) for 
the case MgusY ~ m t- before, we consider the representative choices a^(Mz) — 0.11 and 
as(Mz) = 0.12. For each a 3 (M^) choice, we consider a grid of tan/3 values, holding |V^,(t?^)| 
and m c fixed. We then choose input values for m t and nib (given a^Mz), tan/3, \V c b\, m c ) 
in terms of which all running parameters are uniquely specified at m t : X t (m t ), Xb{m t ) and 
X T {m t ) are given by Eq. (|34]), oti{m t ) are determined by Eqs. (7) and (||) using the central 
values in Eq. (6) R c / t is given by Eq. (|70]), and \V c b\ at scale m t is an input. After integrating 
the RGEs from m t to Mq we check the constraints 

X b (M G ) = X T (M G ) , (75a) 
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^R c/t {M G ) = \V cb {M G )\. (75b) 

If either of these conditions is not satisfied to within 0.2%, we choose another input value 
for mt and mb and repeat the integration. 

We also carry out the RGE calculations with a SUSY scale at 1 TeV. This is done 
exactly as described in the previous section. In addition to the other parameters, we choose 
an input value for the quartic Higgs coupling A at scale m t . We then integrate the two- 



loop standard model RGEs to the SUSY scale and require that Eq. (f|I|) hold to within 
0.1%. For such solutions we apply the other appropriate boundary conditions (given by 
Eqs. (|46|) -(f49D) and integrate the two-loop SUSY RGEs to the GUT scale, where we require 



that X b (M G ) = X T (M G ) and ^jR c/t (M G ) = \V ch {M G )\ to within 0.2%. In our calculation we 
require that rrib,m c and \V c b\ be within the experimentally determined 90% confidence levels 
of the quark mass determinations of GL (4.1 < m& < 4.4 GeV, 1.19 < m c < 1.35 GeV) and 
the recent Particle Data Book value JTTJ] for |V c6 | (0.032 < |V c6 | < 0.054). 



In Fig. 9 the contours of constant \V c b\ are shown in the mt, tan (3 plane for a fixed 
m c = 1.27 GeV. In Figs. 10 and 11 we show the contours obtained by applying only 
the constraint in Eq. (75a) as in Fig. 5 along with the contours obtained by applying both 
Eqs. (75a) and (75b) for fixed m c as in Fig. 9. In Fig. 10 the value of m c is fixed at 1.27 GeV 
and contours of | V c b\ are shown. In Fig. 11 \ V c b\ is fixed at its maximum allowed experimental 
value of \V c b\ = 0.054 (at 90% C.L.) and three values of m c are plotted (corresponding to the 
central m c value and the 90% C.L. values from GL). 

For large tan/3 the effects of including Xb and A r in the RGEs increase \V c b\. In order 
to satisfy \V c b\ < 0.054, the maximum allowed value of tan/? for a 3 (M^) = 0.11 is about 
50(60) for Msusy = m 4 (lTeV); see Fig. 11. For this value of a 3 {M z ) the HRR/DHR model 
predicts that \V c b\ still lies at the upper end of its allowed 90% confidence level range when 
the effects of A& and A T at large tan/3 are included in the two-loop RGEs; see Fig. 10. 
Allowing mfc to become larger than the narrow window mt, = 4.1 — 4.4 GeV requires bigger 
\V c b\ which is unacceptable. The higher b mass contour m b = 5 GeV is not consistent with 
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the GUT scale ansatz for a$(Mz) = 0.11. The largest consistent values of m;, are given in 
Table 4. 



MsTJSY 


a 3 (M z ) 


0.11 


0.12 


m t 


4.56 


5.28 


1 TeV 


4.70 


5.33 



Table 4: Maximum values of 7775 (m&) in GeV consistent with 
the 90% confidence levels of \V c b\ and m c {m c ). 



With a^iMz) = 0.12, \V c b\ can be much closer to its central value, enhancing the plausi- 
bility of the HRR/DHR model, with the only caveat being that low nib (^4.2 GeV) values 
produce \ t {Mo) values which are close to being non-perturbative for most values of tan/?: 
see Figs. 6b, 6d. Notice that the dominant effect of taking the larger value of as(Mz) indi- 
cated by two-loop evolution is to increase the QCD-QED scaling factor i] c , thereby allowing 
\V c b\ to be smaller and in better agreement with experiment. 

Imposing the constraints on m&, m c and \V c b\ also gives the lower limits on the top quark 
mass since the |V^J contours in the smaller tan/3 region are steeper and eventually cross the 



rrib/m T contours ||18||. These lower limits on mt are summarized in Table 5. 



M-SUSY 


a 3 (M z ) 


0.11 


0.12 


m t 


155 (1.45) 


118 (0.75) 


1 TeV 


151 (1.16) 


116 (0.64) 



Table 5: Minimum values of mt(m t ) (tan/?) in GeV consistent with 
the 90% confidence levels of m^m^), \V c b\ and m c (m c ). 



The constraints on rrib/m T , \V c b\ and m c completely determine the allowed region in the 
m t ,tan/3 plane of the HRR/DHR model. Other constraints such as the e parameter for CP 
violation in the neutral kaon system, B mixing or the lighter quark masses affect only the 



other parameters in the model [18 



If the Yukawa unification is assumed to occur at a scale higher than the gauge couplings, 
then the predicted value for \V c b\ will be lower [|J and easier to reconcile with the experimental 
data. 

B. The Giudice Model 

Giudice has proposed a different Yukawa mass ansatz of the form 
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This model uses a geometric mean relation m c = m u m t at the GUT scale to eliminate one 
parameter in the up quark Yukawa matrix. The down quark Yukawa matrix must then 
generate the mixing between the second and third generations to get a value for \V c b\ that 
agrees with experiment. Giudice sets the parameter n in the above mass matrices to be two. 
We see no a priori reason to suppose that this parameter must be an integer and treat it as 
a free parameter. 

We find the generalized one-loop solutions (neglecting and A T in the RGEs) 



IK. 



<m P 



I 25 m e 4n 2 m M r\ T " 
2 m M 9 m T rj fl/ 



|VJ = f=S* 1- = + 
y 2 m c 

Vub\ — , 



(n 2 - 3) rj T 



9 m T 7) fl/ 



(78) 
(79) 
(80) 



28 



Vc m t 

m d = 3m e 1-8 h— -— , (82) 

x rj e \ m M 9 m T %J 

m s = Ollk^hL 1 1 + 8 ^ - ^.TJL±) , (83) 
x rjfj, 3 \ m M 9 m r ?7 M y 

= |/ m r . (84) 

x r\ r 

Notice that at one-loop level to leading order in the mass heirarchy the running \V c b\ is 
related to the strange and bottom Yukawa couplings by 

\VM\ = nR. /b (ji) = n^f\. (85) 



Eqs. (jfo|)-(p^) can be compared to Eqs. (p8|)-(|63|) for the HRR/DHR model, except that 
we have retained the highest non-leading order corrections only for the Giudice model. When 
n = 2 the predicted value of \V c b\ agrees well with the experimental value. On the other 
hand |V^ S | is just at the lower limit of its 90% confidence level. The overall situation can be 
improved somewhat by allowing n to be slightly larger than two. 

The leading term in Eq. ( [78] ) can be recognized as the Oakes relation jH| between the 



Cabibbo angle and the quark masses, tan# c ~ ym d /m S} supplemented by the Yukawa 
unification relation m d jm s = Qrrie/m^. Notice that this relation involving the first and 
second generations does not run, so the prediction of the Cabibbo angle is insensitive to the 
size of the gauge and Yukawa couplings. The two-loop effects for the most part increase 
«3 and hence the QCD scaling factors r\ q . The influence of two-loop contributions in the 
running of the Yukawas is small. 

For tan/3^10, and A r can be neglected in the RGEs; then the relation for m u in 
Eq. ([81]) implies an upper limit on m t |7J . However further solutions for jm T are possible 
with large tan/?, as can be seen in Figure 5. In the allowed m\ ) jm T band at large tan/? the 
predicted value for m u from Eq. ( |S~1"D is still satisfactory, since mt is in the same range as 
found for the small tan (3 solutions. 

The CP-violating phase is not very well constrained in the Giudice model since the 
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phase does not enter in the well-measured CKM elements; in fact the phase can assume 
almost any non-zero value within its zero to 2n range. Correspondingly CP asymmetries 
to be measured in B decays are not very constrained in the model ||44|| . In contrast, the 
CP- violating phase in the HRR/DHR model is almost uniquely determined by \V US \ and the 
CP- violating asymmetries are predicted precisely. This remain the case at the two-loop level. 
In the HRR/DHR scheme the dependence on a^iMz) cancels out in quark mass ratios, and 
since the constraint on the phase arises from the first and second generation mixing angles, 
there is no dependence of the phase on \ t . 

V. CONCLUSION 

We have investigated unification scenarios in supersymmetric grand unified theories using 
the two-loop renormalization group equations. Our primary conclusions are the following: 

(1) Given the experimentally determined values for oti and a 2 at M z , the RGEs pre- 
dict 0:3 ~ 0.111(0.122) at one-loop (two-loop) for M SUSY = m t and a 3 ~ 0.106(0.116) for 
M SUSY = 1 TeV. Including the Yukawa couplings in the two-loop evolution of the gauge cou- 
plings decreases a 3 (M z ) by only a few percent. Thus the values of a 3 (M z ) — 0.12 obtained 
experimentally at LEP II are also theoretically preferred if GUT scale thresholds effects or 
intermediate scales are not important. 

(2) For any fixed value of a 3 (M z ) and rrib there are just two allowed solutions for tan/? 
for a given top mass if m t ^180 GeV; the larger solution has tan/5 > m t /m b and the smaller 
solution is sin/5 ~ 0.78(m i /150GeV). Allowing for some uncertainty in a 3 (M z ), and 
M SUSY , these unique solutions for tan/3 at given m t become a narrow range of values. For 
m t ~ 180 — 200 GeV the value of tan/5 changes rapidly with m t . 

(3) With At,, A T unification we find an upper limit mt^200 GeV on the top quark mass by 
requiring the successful prediction of the m\ ) jm T ratio; we also obtain lower limits m t ^150 
GeV (115 GeV) for a 3 (M z ) = 0.11(0.12) from evolution constraints on m b , m c and |V^,|. 
These lower limits are only mildly sensitive to M SUS y- 
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(4) The effects of raising M SUSY is to decrease both a G and M G and to decrease the values 
of a 3 (M z ) that yields successful unification. Also the allowed band for the m h /m T ratio in 
the m t , tan (5 plane is shifted towards slightly higher top masses. This in turn slightly reduces 
the prediction for (V^l in models that utilize the relation ^ \ C (M G ) / \ t (M G ) = \V c b(M G )\. 

(5) In the HRR/DHR model we find an upper limit on the supersymmetry parameter 
tan ,3^50(60) for M SUSY = m t (lTeV) if a 3 {M z ) ~ 0.11; for a 3 {M z ) = 0.12 the solutions at 
large tan/5 extend into the region for which \b(M G ) is non-perturbative. 

(6) For the value a 3 (M z ) ^ 0.12 indicated by the two- loop RGEs, the agreement of 
the \V c b\ prediction of the HRR/DHR ansatz with experiment is improved. In fact for 
az(M z ) = 0.12 and M SUSY = 1 TeV the central values for jV^I and the mass ratio m b /m T 
almost coincide in the m t ,tan(3 plane; see Fig. lOd. This result is more general than the 
HRR/DHR ansatz, and occurs for any model with the GUT scale relation \V c b\ = \J A c / A t . 

(7) With at(M z ) ~ 0.12 a large top Yukawa coupling is needed to achieve the correct 
rrib/m T ratio, and the theory is in some jeopardy of having a non-perturbative A t (M G ) if m b 
is smaller than about 4.2 GeV. 

(8) GUT unification of A T , \ b and \ t can be realized for tan/3^50. 

(9) The predictions for the CP asymmetries in the HRR/DHR model are largely unaf- 
fected by our two-loop analysis. 

(10) We have found new solutions to the Giudice model for large tan/5. These results 
require the inclusion of \b and A T in the RGEs, and therefore could not be obtained in 
Giudice's analytic treatment at one-loop. 
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VI. APPENDIX 



To consider a specific ansatz for Yukawa matrices at the GUT scale at the two-loop level 
requires knowledge of the RGEs. These can be derived from formal expressions that exist 
in the literature For the supersymmetric model with two Higgs doublets, the one- [2~3 
and two-loop RGEs can be written for general Yukawa matrices as 



dgj 
dt 
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These equations agree with those in the last paper in Ref. p5[ for the case where the Yukawa 
matrices are diagonal, if the following minor corrections are made: (1) 631 should be decreased 
by a factor three; (2) the parenthesis in the second term of 7^ should come before the a 2 ; 
(3) the first term of should have a factor a\ instead of a 2 .. 
The two-loop RGEs for the standard model are [E6f 
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(111) 



These renormalization group equations are those given in the classic papers of Machacek and 
Vaughn after replacing H — ► U^, F D — > D^, F L — > E^, and making the following corrections 
to Eq. ( |101| ) mentioned in the paper of Ford, Jack and Jones PB| : (1) The \g% term in the 
one-loop beta function has a coefficient 9 instead of 1. (2) The \g\g\ term in the two- loop 
beta function has a coefficient +117/20 instead of —117/20. (3) The Xgf in the two-loop 
beta function has a coefficient +1887/200 instead of —1119/200. 



36 



REFERENCES 



[1] For review and references see e.g. R. Mohapatra, Unification and Super symmetry, 2nd 
Edition, Springer- Verlag (1992); G. G. Ross, Grand Unified Theories, Addison- Wesley 
(1984). 

[2] U. Amaldi, Wim de Boer, and H. Furstenau, Phys. Lett. B260, 447 (1991); J. Ellis, 
S. Kelley and D. V. Nanopoulos, Phys. Lett. B260 131 (1991); P. Langacker and M. Luo, 
Phys. Rev. D44, 817 (1991). 

[3] U. Amaldi, et al., Phys. Rev. D36, 1385 (1987); G. Costa, J. Ellis, G. L. Fogli, 
D. V. Nanopoulos and F. Zwirner, Nucl. Phys. B297, 244 (1988); J. Ellis, S. Kelley and 
D. V. Nanopoulos, Phys. Lett. B249, 442 (1990); S. Kelley, J. L. Lopez, D. V. Nanopou- 
los, H. Pois and K. Yuan, Phys. Lett. B273, 423 (1991); F. Anselmo, L. Cifarelli, A. Pe- 
termann and A. Zichichi, CERN preprint CERN-PPE-9 1-123 (1991); J. Ellis, S. Kelley 
and D. V. Nanopoulos, Phys. Lett. B287, 95 (1992). 

[4] H. Arason, D. J. Castano, B. Keszthelyi, S. Mikaelian, E. J. Piard, P. Ramond and 
B. D. Wright, Phys. Rev. Lett. 67, 2933 (1991), P. Ramond, University of Florida 
preprint UFIFT-92-4; H. Arason, D. J. Castano, E. J. Piard and P. Ramond, University 
of Florida preprint UFIFT-92-8. 

[5] H. Arason, D. J. Castano, B. Keszthelyi, S. Mikaelian, E. J. Piard, P. Ramond and 

B. D. Wright, University of Florida preprint UFIFT-91-33 (1991). 

[6] S. Dimopoulos, L. J. Hall and S. Raby, Phys. Rev. Lett. 68, 1984 (1992); Phys. Rev. 
D45, 4192 (1992). 

[7] G. F. Giudice, University of Texas preprint UTTG-5-92 (1992). 

[8] B. R. Greene, K. H. Kirklin, P. J. Miron and G. G. Ross, Nucl. Phys. B292, 606 (1987); 
J. L. Lopez and D. V. Nanopoulos, Phys. Lett. B251, 73 (1990); J. A. Casas and 

C. Munoz, Nucl. Phys. B322, 189 (1990), F. del Aguila, M. Masip and L. A. C. P. da 

37 



Mota, CERN preprint CERN-TH.6233/91 (1991); K. S. Babu and Q. Shan, Bartol 
Research Institute preprint BA-92-27 (1992). 

[9] G. G. Ross and R. G. Roberts, Nucl. Phys. B377, 571 (1992); R. Arnowitt and R Nath, 
Phys. Rev. Lett. 69, 725 (1992); S. Kelley, J. L. Lopez, D. V. Nanopoulos, H. Pois and 
K. Yuan, CERN preprint CERN-TH.6498/92 (1992). 

[10] See e.g. S. Bethke and S. Catani, CERN-TH.6484/92 (1992); S. Bethke, talk at XXVI 
Int. Conf. on HEP, Dallas (1992). 

[11] Particle Data Book, Phys. Rev. D45 (1992). 

[12] J. Ellis, D. V. Nanopoulos and S. Rudaz, Nucl. Phys. B202, 43 (1982); R. Arnowitt 
and P. Nath, Northeastern University preprint NUB-TH-3051-92 (1992); J. L. Lopez, 
D. V. Nanopoulos and A. Zichichi, CERN preprint CERN-TH.6554/92 (1992); J. Hisano, 
H. Murayama and T. Yanagida, Phys. Rev. Lett. 69, 1014 (1992); Tohoku University 
preprint TU-400 (1992). 

[13] H. Georgi and S. Glashow, Phys. Rev. Lett. 32, 438 (1974); H. Georgi, H. Quinn and 
S. Weinberg, Phys. Rev. Lett. 33, 451 (1974); S. Dimopoulos and H. Georgi, Nucl. Phys. 
B293, 150 (1981); N. Sakai, Z. Phys. Cll, 153 (1981). 

[14] M. Chanowitz, J. Ellis and M. Gaillard, Nucl. Phys. 128, 506 (1977); A. Buras, J. Ellis, 
M. Gaillard and D. Nanopoulos, Nucl. Phys. 135, 66 (1978). 

[15] M. B. Einhorn and D. R. T. Jones, Nucl. Phys. 196, 475 (1982); J. Ellis, D. V. Nanopou- 
los and S. Rudaz, Nucl. Phys. 202, 43 (1982). 

[16] J. Harvey, P. Ramond and D. B. Reiss, Phys. Lett. B92, 309 (1980); Nucl. Phys. B199, 
223 (1982). 

[17] S. Kelley, J. L. Lopez and D. V. Nanopoulos, Phys. Lett. B274, 387 (1992); J. Ellis, 
S. Kelley, D. V. Nanopoulos, Nucl. Phys. B373, 55 (1992). 



38 



[18] V. Barger, M. S. Berger, T. Han and M. Zralek, Phys. Rev. Lett. 68, 3394 (1992). 
[19] R. J. Oakes, Phys. Lett. B29, 683 (1969). 

[20] H. Fritzsch, Phys. Lett. B70, 436 (1977); Phys. Lett. B73, 317 (1978); Nucl. Phys. 
B155, 189 (1979). 

[21] H. Georgi and C. Jarlskog, Phys. Lett. B86, 297 (1979). 

[22] B. Stech, Phys. Lett. B130, 189 (1983); M. Shin, Phys. Lett. B145 285 (1984); 
M. Gronau, R. Johnson and J. Schechter, Phys. Rev. Lett. 54, 2176 (1985); H. Fritzsch, 
Phys. Lett. B166, 423 (1986); M. Soldate, M. H. Reno and C. T. Hill, Phys. Lett. 
B179, 95 (1986); R. Johnson, J. Schechter and M. Gronau,, Phys. Rev. D33, 2641 
(1987); M. Tanimoto and A. Takatsu, Phys. Rev. D34, 1474 (1986); F. del Aguila, 
G. L. Kane and M. Quiros, Phys. Lett. B196, 531 (1987); C. H. Albright, C. Jarl- 
skog and B. Lindholm, Phys. Lett. B199, 553 (1987); Phys. Rev. D38, 872 (1988); 
M. Gronau, R. Johnson and J. Schechter, Phys. Lett. B201, 151 (1988); C. H. Albright, 
Phys. Lett. B224, 213 (1989); H. P. Nilles, M. Olechowski and S. Pokorski, Phys. Lett. 
B228, 406 (1989); X. G. He and W. S. Hou, Phys. Rev. D41, 1517 (1990); P. Kaus 
and S. Meshkov, Phys. Rev. D42, 1863 (1990); C. H. Albright, Phys. Lett. B246, 451 
(1990); D. Ng and Y. J. Ng, Mod. Phys. A6, 2243 (1991); J. L. Rosner and M. Worah, 
Enrico Fermi Institute preprint EFI92-12 (1992). 

[23] B. Pendleton and G. G. Ross, Phys. Lett. B98, 291 (1981); N. Cabibbo, L. Maiani, 

G. Parisi and R. Petronzio, Nucl. Phys. B158, 295 (1979); C. T. Hill, Phys. Rev. D24, 
691 (1981); E. A. Paschos, Z. Phys. C26, 235 (1984); J. W. Halley, E. A. Paschos and 

H. Usler, Phys. Lett. B155, 107 (1985); J. Bagger, S. Dimopoulos and E. Masso, Nucl. 
Phys. B253, 397 (1985); M. Tanimoto, T. Hayashi, R. Najima and S. Wakaizumi, Prog. 
Theor. Phys. 76, 1098 (1986); K. S. Babu and E. Ma, Europhysics Lett. 3, 437 (1987); 
M. Tanimoto, Y. Suetake and K. Seuba, Phys. Rev. D36, 2119 (1987); C. H. Albright 
and M. Lindner, Phys. Lett. B213, 347 (1988); Z. Phys. C44, 673 (1989); E. M. Freire, 



39 



G. Lazarides and Q. Shan, Mod. Phys. Lett. A5, 2453 (1990); G. Cvetic and C. S. Kim, 
Dortmund University preprint DO-TH-92-08 (1992). 

[24] K. Inoue, A. Kakuto, H. Komatsu and S. Takeshita, Prog. Theor. Phys. 67, 1889 (1982). 

[25] D. R. T. Jones and L. Mezincescu, Phys. Lett. B136, 242 (1984); P. West, Phys. Lett. 
B137, 371 (1984); A. Parkes and P. West, Phys. Lett. B138, 99 (1984); D. R. T. Jones 
and L. Mezincescu, Phys. Lett. B138, 293 (1984); J. E. Bjorkman and D. R. T. Jones, 
Nucl. Phys. B259, 533 (1985). 

[26] M. E. Machacek and M. T. Vaughn, Nucl. Phys. B222, 83 (1983); Nucl. Phys. B236, 
221 (1984); Nucl. Phys. B249, 70 (1985); M. Fischler and J. Oliensis, Phys. Lett. B119, 
385 (1982); Phys. Rev D28, 2027 (1983); I. Jack, J. Phys. A16, 1083 (1983); I. Jack and 

H. Osborn, J. Phys. A16, 1101 (1983); I. Jack, Phys. Lett. B147, 405 (1984); I. Jack 
and H. Osborn, Nucl. Phys. B249, 472 (1985); C. Ford, I. Jack and D. R. T. Jones, 
NSF-ITP-92-21 (1992). 

[27] P. Langacker, University of Pennsylvania preprint UPR-0492T (1992). 

[28] A. Sirlin, Nucl. Phys. B332, 20 (1990); P. Langacker, Phys. Rep. 72, 185 (1981); 
A. Buras, Rev. Mod. Phys. 52, 199 (1980); G. 't Hooft, Nucl. Phys. B61, 455 (1973). 

[29] W. Marciano, Phys. Rev. D29, 580 (1984). 

[30] A. Giveon, L. J. Hall and U. Sarid, Phys. Lett. B271, 138 (1991). 

[31] H. E. Haber and R. Hempfling, Santa Cruz preprint SCIPP-91/33 (1992). 

[32] S. Weinberg, Phys. Lett. B91, 51 (1980). 

[33] L. J. Hall, Nucl. Phys. B178, 75 (1981). 

[34] R. Barbieri and L. J. Hall, Phys. Rev. Lett. 68, 752 (1992). 

[35] K. Hagiwara and Y. Yamada, KEK-TH-331 (1992); A. E. Faraggi, B. Grinstein and 



40 



S. Meshkov, SSCL-Preprint-126 (1992). 

[36] O. V. Tarasov, A. A. Vladimirov and A. Y. Zharkov, Phys. Lett. B93, 429 (1980); 
S. G. Gorishny, A. L. Kateav and S. A. Larin, Sov. J. Nucl. Phys. 40, 329 (1984); 
S. G. Gorishny et al, Mod. Phys. Lett. A5, 2703 (1990). 

[37] R. Tarrach, Nucl. Phys. B183, 384 (1981). 

[38] In the notation of Ref. || these parameters are y - 1 = (A t (M G )/A t ) 1/f y /12 . In the nota- 
tion of Ref. 0, y = £ and x = (7£>/7e) _1 V . 

[39] J. Gasser and H. Leutwyler, Phys. Rep. 87, 77 (1982). Referred to as GL in the text. 

[40] E. Ma and S. Pakvasa, Phys. Lett. B86, 43 (1979), Phys. Rev. D20, 2899 (1979); 
K. Sasaki, Z. Phys. C32, 149 (1986); B. Grzadkowski and M. Lindner, Phys. Lett. 
B193, 71 (1987); M. Olechowski and S. Pokorski, Phys. Lett. B257, 388 (1991). 

[41] K. S. Babu, Z. Phys. C35, 69 (1987); B. Grzadkowski and M. Lindner, Phys. Lett. 
B198, 64 (1987); See also the last paper of Ref. |0 . 



[42] CERN library integration routine DDEQMR based on the algorithm in N. Lance, Nu- 
merical Methods for High Speed Computers, Iliffe and Sons, London (1960), 56. 

[43] V. Barger, M. S. Berger, A. L. Stange and R. J. N. Phillips, Phys. Rev. D45, 4128 
(1992); V. Barger, K. Cheung, A. L. Stange and R. J. N. Phillips, Madison preprints 
MAD/PH/696, MAD/PH/704 (1992); H. Baer, M. Bisset, C. Kao and X. Tata, Phys. 
Rev. D46, 1067 (1992); J. F. Gunion, R. Bork, H. E. Haber and A. Seiden, Santa Cruz 
preprint SCIPP-91-34 (1991); J. F. Gunion, H. E. Haber and C. Kao, University of 
California-Davis preprint UCD-91-32 (1991); J. F. Gunion and L. H. Orr, University 
of California-Davis preprint UCD-91-10 (1991); H. E. Haber, R. Hempfhng and Y. Nir, 
Santa Cruz preprint SCIPP-91-40 (1992); Z. Kunszt and F. Zwirner, ETH preprint 
ETH-TH-91-7 (1991). 



41 



[44] S. Dimopoulos, L. J. Hall and S. Raby, Lawrence Berkeley Lab preprint LBL-31431 
(1992); U. Sarid and Y. Nir, Lawrence Berkeley Lab preprint LBL-32563, (1992). 



42 



Figures 



Fig. 1. Allowed GUT parameter space for m t = 150 GeV with (a) M SUSY = mt (one- 
loop RGE) (b) M SUSY = m t (two-loop RGE) (c) M SUSY = 1 TeV (one-loop RGE) (d) 
M-susy = 1 TeV (two-loop RGE) versus the running mass scale /i. The shaded region 
denotes the range of GUT coupling and mass consistent with the la ranges of a\[Mz) and 
a<i(Mz)\ the curves for a^/u,) represent extrapolations from the GUT parameters. We have 
omitted the contributions from Yukawa effects here which depend on tan j3. 

Fig. 2. Gauge coupling unification with two-loop evolution for (a) M SUSY = m t (b) M SUSY = 
1 TeV taking m t = 150 GeV and neglecting Yukawa couplings; fi is the running mass scale. 

Fig. 3. The QCD-QED scaling factors r/f of Eq. fllTD are shown for f = s,c,b versus as(Mz), 
assuming running quark masses mjimj) of m t = 170 GeV, m h = 4.25 GeV, m c = 1.27 GeV. 

Fig. 4. The top Yukawa coupling at the GUT scale determined at the one-loop level is 
plotted versus a 3 (M z ) for m t = 170 GeV and m& = 4.25 GeV. 

Fig. 5. Contours of constant m& in the mt, taxi /3 plane obtained from the RGEs with (a) 
Msusy = mt, a 3 (M z ) = 0.11; (b) M SUSY = m t , a 3 (M z ) = 0.12; (c) M SUSY = 1 TeV, 
as(Mz) = 0.11; (d) M SUSY = 1 TeV, a^iMz) = 0.12. The shaded band corresponds to 
the 90% confidence level range of m;, from Ref. |39| (m& = 4.1 — 4.4 GeV); the dotted curve 
corresponds to m& = 5.0 GeV. The curves shift to higher m t values for increasing a 3 (M^) or 
increasing M SUSY . 
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Fig. 6. The Yukawa couplings X^Mg) and A&(Mg) = X t (Mg) at the GUT scale with (a) 
m susy = m u a 3 {M z ) = 0.11; (b) M SUSY = m t , a 3 {M z ) = 0.12; (c) M SUSY = 1 TeV, 
a 3 (M z ) = 0.11; (d) M SUSY = 1 TeV, a 3 (M z ) = 0.12. The Yukawa couplings become larger 
for higher a 3 (M z ) or higher M SUSY . The perturbative condition A^3.3 from Eq. ([52]) is 
satisfied except for the lowest b mass value m b = 4.1 GeV for a 3 {M z ) = 0.12. The solid dots 
denote A r = \ b = At unification. 

Fig. 7. Two- loop evolution of the Yukawa couplings (a) At(/i) (b) \ b (fi), A r (/x) from low 
energies to the GUT scale for the case a 3 (M z ) = 0.12 and M SUSY = 1 TeV. We take tan f3 = 
20 and the values of m t = 198, 197, 196, 181 GeV specified by the m b = 4.1, 4.25, 4.4, 5.0 
GeV contours in Fig 5d. 

Fig. 8. Two-loop evolution of the quark Yukawa ratio R c / t = A c /Af and the CKM matrix 
element \V c b\ for (a) M SUSY = m t and (b) M SUSY = 1 TeV. We have taken «3 = 0.11, tan (3 = 
5 and have chosen the top and bottom quark masses such that \J R c /t(Mc) = \ V c },(Mg)\ and 
m c = 1.27 GeV: (a) \V cb (m t ) \ = 0.054, m t = 180 GeV, m b = 4.33 GeV; (b) \V cb (m t ) \ = 0.050, 
m t = 189 GeV, m b = AAA GeV. 

Fig. 9. Contours for constant \V cb \ at fixed m c = 1.27 GeV in the mt, tan /3 plane obtained 
from the RGEs with (a) M SUSY = m t , a 3 {M z ) = 0.11; (b) M SUSY = 1 TeV, a 3 {M z ) = 0.12. 

Fig. 10. Comparison of contours for constant \V cb \ and constant m b in the mt, tan /3 plane 
from the RGEs, taking m c = 1.27 GeV, for (a) M SUSY = m t , a 3 (M z ) = 0.11; (b) M SUSY = 
m t , a 3 [M z ) = 0.12; (c) M SUSY = 1 TeV, a 3 {M z ) = 0.11; (d) M SUSY = 1 TeV, a 3 [M z ) = 
0.12. The shaded band indicates the region where the 90% confidence limit is satisfied for 
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m&. The right-most contours are discontinued when X t (Mo) exceeds 6. 

Fig. 11. Comparison of contours for constant m c and constant m h in the m t , tan j3 plane from 
the RGEs, taking \V c b\ equal to its upper limit 0.54, for (a) M SUSY = m t , a^{Mz) = 0.11; 
(b) M SUSY = mt, a 3 (M z ) = 0.12; (c) M SUSY = 1 TeV, a 3 (M z ) = 0.11; (d) M SUSY = 1 
TeV, a 3 (M z ) = 0.12. The shaded band indicates the region where 90% confidence limits 
are satisfied for all three constraints: m b , m c and An X marks the lower limit of this 

shaded band and corresponds to the values in Table 5. 
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